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ABSTRACT 



Exact analytical formulae are derived for the potential and mass ratio as a 
function of Lagrangian points position, in the classical Roche model of the close 
binary stars. 



Subject headings: stars: rotation — binaries: close 



1. 



Introduction 



The Roche model is widely used in interpretation of the close binary star observations. 
Several authors derived the various approximations to solve the Roche problem and presented 
numerical tables, see e.g. Plavec & Kratochvil (1964); Kippenhahn & Thomas (1970); 
Eggleton (1983); Mochnacki (1984); Morris (1985, 1994). Not intending to solve the whole 
problem analytically, I rather show here the way how some analytical relations can be found. 
The idea is to reverse the problem: instead of finding, e.g., the first Lagrangian point, xi, as 
function of binary mass ratio q we seek solution for q as function of x\. 



The basic equation of the classic Roche problem is the formula for surfaces of the 
(primary) star as equipotential : 



2. Basic equation 



^{x,y,z) = ( x - j-p- J +y 2 + 




(1 + q) (-1 + xf + y 2 + z 2 



2q 
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(1 + q) \J 1 x 2 + y 2 + z 2 
(1) 
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I use in eq. (1) notations of Mochnacki (1984): the x axis is aligned along the stars' centers, z 
axis is parallel to the rotation axis, \& is "normalized" potential, in units of G(M 1 +M 2 )/(2 A), 
q = M 2 /Mi < 1, and x,y,z are in units of A, the distance between the centers of binary 
components with masses Mi and M 2 . As I here consider only the first (inner) and second 
(outer) Lagrangian points problem, it is sufficient to consider the eq. (1) at the x axis. 



3. The first Lagrangian point 

We write down eq. (1) with y = z = 0, < x < 1, and look for the minimum of the 
function: 

Hf^x) = % + - 2 +( x - -±-) \ (2) 

W (l + q)(l-x) (l + q)x V 1 + 9/ 

at some x — X\, with < Xi < 1. The important observation, from eq. (2), is that 
^i(q,x) = ^1(1/5,1 — x), if < x < 1 (not in general case!). From eq. (2), we find the 
condition d^i(x)/dx = 0, which we rewrite after some algebra as function of q(xi): 

q{Xl) = xf(3-3x 1 + xj) ■ (3) 

We notice, from eq. (3), the elegant relation (also having a clear physical meaning) q(xi) q(l — 
x\) = 1. We underline that eq. (3) gives the fully analytical and exact relation between values 
of q and of the first Lagrangian point x\. See, e.g. table 1 in (Mochnacki 1984), where Q 
and XI stand for our q and xi, respectively, and compare the difference in difficulty of 
calculations by method of Mochnacki (1984) and by formulas (3,4,6,9) of this note. 

Now in order to find value of the potential corresponding to the first Lagrangian point, 
we may use the eqs. (2) and (3) together, or, in the spirit of this note, exclude q from eqs. 
(2) and (3) and find explicit function \& 1 (a; 1 ): 

T . . 3- 12t+ 15t 2 - 10t 3 - 4t 4 

Vi(xi) = 5 :t = xi(l-x 1 ). (4) 

U ; (_l+2t + t2) 2 1 ; 

We introduced the additional variable t into Eq. (4) in order to explicitly show that ^i(xi) = 
^1(1 — xi) if < x\ < 1 (not in general case!). Eq. (4) gives the fully analytical and exact 
relation between the values of the first Lagrangian point x\ and of the corresponding potential 
^>i(xi). See, e.g. table 1 in (Mochnacki 1984), where CI and XI stand for our ^ r i(xi) and 
x±, respectively. 
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4. The second Lagrangian point 

Now we look for minimum of the function (note the difference from eq. [2]): 

*2(a0 = t 1-7 r + t —s 1- (x — 1 , (5) 

2y ! (l + q)(x-l) (l + q)x V 1 + 9/ 

at some x = x 2 , with 1 < x 2 < 2. Repeating the procedure of the section 3, we get the 
solution for q as function of x 2 . 

«m = { 2 x r 1) \ i ] +X2+x2 j - (6) 

yv 2J x\{2-x 2 ) (l-x 2 + x 2 2 ) y 1 

We notice that by contrast with eq. (3), it is not evident from eq. (6) that we have the 
relation q(x 2 ) q(l — x 2 ) = 1, because here x 2 > 1, but (1 — x 2 ) < 0, while the eq. (6) is 
derived under condition 1 < x 2 < 2; in the case of eq. (3), both X\ and (1 — x±) were in the 
(open) interval (0,1)). 

To prove the validity of the relation q(x 2 ) q{l — x 2 ) = 1, we should return to the basic 
equation (1), put there y = z = 0, x < 0, and look for the minimum of the function: 

*'W - (i + ,hi-x) - (TT^ + { X -T^) ■ (7) 

at some x = X3, with < 0. Note the differences between functions , I' 1 (a;), \l/ 2 (x), and 
^3(2;). Repeating above procedure for ^(x), we get: 

n(rr \ (2 - ar 3 ) x§ (1 - rr 3 + arj) , , 

(x 3 - 1) (l + x 3 + ^) 

Now, from eqs. (6) and (8), we have q(x 2 ) q(xs) = 1, if x 3 = 1 — x 2 , and if x 2 > 1, £3 = 
1 — x 2 < 0, QED. Eq. (6) (together with eq. [8]) gives the fully analytical and exact relation 
between q and the value of the second Lagrangian point x 2 . See, e.g. table 1 in (Mochnacki 
1984), where Q and X2 stand for our q and x 2 , respectively. 

Now in order to find value of the potential corresponding to the second Lagrangian 
point, we may use the eqs. (5) and (6) together, or, again in the spirit of this note, exclude 
q from eqs. (5) and (6) and find explicit function ^2(^2): 

T , , -1 -4x 2 + 27a:| -36xl + 9xt + 18x1 - Ux 6 2 + Axl 

ty 2 (x 2 ) = ^5 -. (9) 

V ; (-l + 2x 2 + x 2 2 -2xl + x 4 2 ) 2 y ' 

Eq. (9) gives the fully analytical and exact relation between the values of the second La- 
grangian point x 2 and of the corresponding potential ^2(^2) • See, e.g. table 1 in (Mochnacki 
1984), where C2 and X2 stand for our ^2(^2) and x 2) respectively. 
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5. Summary 

In this short note we present the exact analytical relations for the first and second 
Lagrangian point of the classical Roche problem. In practice, the "exact" formulas are not 
necessarily the most convenient ones, hence different approximate expressions in literature 
(see (Kopal 1959; Plavec & Kratochvil 1964; Kippenhahn & Thomas 1970; Eggleton 
1983; Mochnacki 1984; Morris 1985, 1994), among others), and some relevant approximate 
formulas will be given elsewhere. Still, the exact formulas have their own beauty and are 
more relevant as solutions to the classical problems such as the problem by Roche (1847). 

My thanks are due to anonymous referee for encouraging criticism. 
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